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AN  EMPIRICAL  MODEL  FOR  NEAR  MILLIMETER  WAVE 
SNOW  EXTINCTION  AND  BACKSCATTER  (U) 

♦Douglas  R.  Brown.  Ph.P . 

OS  Aray  Atmospheric  Sciences  Laboratory 
White  Sands  Missile  Range,  New  Mexico  88002 


_ A  snow  extinction  and  badtscattfcr  model  for  near  millimeter  waves 

was  developed  to  complete  a  set  of  natural  atmosphere  propagation  models 
for  iuclusion  in  the  Electro-Optics  System  Effects  Library  (EOSAEL) 
(£)*  This  model  is  compact,  simple,  fast,  and  uses  consonly  available 
meteorological  data  as  inputs.  The  accuracy  is  commensurate  with  the 
accuracy  of  those  input  parameters;  where  the  input  parameters  dc  not 
properly  characterise  the  environment,  further  input  specification  is 
available.  For  snow,  the  rain  equivalent  accumulation  rate  is  the  funda¬ 
mental  input,  with  air  temperature  providing  the  distinction  between  dry 
and  wet  snowtypes.  For  millimeter  wave  (MMW),  the  classification  of  the 
snowflake  by  its  ice-tf'  -liquid-water  ratio  is  necessary,  due  to  the  large 
differences  in  the  complex  indices  of  refraction  of  ice  and  water,  (2)% 

The  snow  extinction  model  is  presented  first,  including  a  re-iew  of 
the  data  used  to  generate  it.  The  bsckscatter  ~»del  is  than  discussed, 
followed  by  a  comparison  between  the  model  and  measurements  made  at  the 
Cold  Regions  Research  and  Engineering  Laboratory  (CREEL)  SNOW-ONE  teat  in 
Janus ry-February  1982,(3). 

* 

l,  SNOW  EXTINCTION 

While  the  literature  on  the  interaction  of  MMW  with  snow  is  limited, 
the  complete  literature  will  not  be  revietvsd  here;  a  general  review  has 
been  carried  out  by  Kobayaahi  (4).  Only  those  results  used  In  the  model 
will  be  discussed. 

The  general  classification  of  snowtypes  is  discussed  in  detcil  by 
Nlshitsuji  (5).  Four  types  are  described:  dry,  moist,  wet,  and  watery, 
distinguished  by  the  deesity  of  the  snowflake.  Because  no  other  determi¬ 
nation  of  saowflake  density  was  reported  in  the  literature,  it  is  diffi¬ 
cult  to  follow  this  classification  procedure;  in  practice,  saowtypes  were 
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Inferred  by  comparison  with  the  Nishitaujl  extinction  results  where  pos- 
sible.  Data  for  extinction  by  wet  snow  at  15  GHz,  moist  snow  at  35  GHz, 
and  watery  snow  at  50  GHz  are  presented.  A  Hie  scattering  type  model  was 
developed  by  Nishitsuji,  but  Is  cumbersome,  and  only  bis  experimental 
results  were  employed. 

Data  at  11,  15,  24,  and  48  GHz  are  presented  by  Oomori  and  Aoyagl 
(6),  although  no  indication  of  the  saowtypes  is  given.  A  unique  confocal 
resonator  was  used  to  measure  snow  attenuations.  Comparison  with  other 
data  seems  to  indicate  that  moist  snow  propagation  was  measured.  The 
data  are  modeled  with  an  aR^  relation,  which  is  derived  using  a  theory  of 
snowfall  statistics  and  in  which  R  is  rain  equivalent  enowrate  and  a  and 
b  are  parameters.  The  a  parameter  data  agrees  well  with  other  data,  but 
at  24  and  48  GHz  there  is  poor  agreement  for  the  b  parameter. 


_ Moist  snow  data  is  provided  at  35  GHz  by  Robinson  (7).  The  data  set 

is  of  limited  extent,  having  only  seven  useable  values.  Kalinkin  et  al 
(8)  measured  ury  snow  at  35  (Hz,  with  but  a  slJ'ht  improvement:  eight 
points. 

Only  two  data  points  exist  in  the  literature  for  frequencies  above 
50  GHz:  one  at  140  GHz  for  vet  snow  by  Richard,  Hammerer,  and  Reitz  (9), 
and  one  at  312.5  CHz  for  dry  snow  by  Babkin  at  al  (10).  The  Richard 
paper  also  reports  visual  range  for  the  same  snow  data.  Readers  should 
cote  the  error  in  reference  (9),  figure  25,  where  the  a  parameter  is 
written  as  1.37  but  clearly  should  be  3.7  Instead. 

Analysis  of  the  data  shows  that  an  &R*5  relation  for  snow  is  the 
simplest  iotu  for  a  scow  model,  but  dcea  such  a  model  have  theoretical 
foundations? 
guidelines. 


A  review  of  aR°  models  for  rain  extinction  provides  the 


Relations  of  the  form  aR^  can  be  derived  for  rain  in  the  Rayleigh 
and  optical  limits.  This  result  relies  on  the  use  of  negative  exponen¬ 
tial  or  modified  gamma  dropslze  distribution  and  a  power  series  represen¬ 
tation  for  the  particle  forward  scattering  amplitude.  Where  such  power 
series  representations  cannot  be  found,  direct  fits  to  Mie  scattering 
calculations  are  used.  Such  models  provide  highly  accurate  results  for 
only  n  limited  range  of  rainrates,  but  application  over  wider  ranges  is 
useful  because  of  the  general  dominance  of  errors  introduced  from  the 
beginning  by  the  arbitrary  raindrop  size  distribution.  Errors  can  ba 
significant,  but  may  be  minimized  by  the  judicious  selection  of  raludrop 
distribution  (11). 
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- Rrrnra  in  the  anew  model  can  be  expected  to  be  even  larger.  Errors 

of  drop  or  snowflake  distribution  are  compounded  by  the  snowflake  orien¬ 
tation  dependence  of  the  forward  scattering  amplitude,  which  now  is  due 
to  a  highly  asymmetric  flake  or  aggregate*  Some  of  this  asymmetry  will 
be  reduced  by  volume  averaging  over  the  falling  snowflake  orientations, 
but  larga  variations  due  to  snowflake  type  oust  be  anticipated*  The 
similarity  or_ lark  thereof -between  the  a  and  b  parameter  frequency  and 
particle  size  dependencies  will  indicate  how  complete  this  averaging 
process  is*  For  snowflakes  that  are  truly  amorphous,  with  no  orientation 
forces  in  frei  fall,  t-Mn  toehn-tgnp  «hmt~M  provide  a  reasonable  model. 
This  cancut  be  expected— to  he  true  for  all  anowtypea*  If  flake  orien¬ 
tations  are  r^t  random,  then  aignifir.ant  fluctuations  in  extinction  and 
backscatter  from  the  model  will  occur.  Because  this  model  neglects 
nelarization,  it  oust  also  be  used  cautiously,  since  polarization  effects 
I  rpn„he  .expended  tn  he  large  for  anw.  The  treatment  of  polarization  and 

— snowflake — orientation  are _ heyand  the  scope  of  a  model  of  this  type. 

[  There  are  probably  few,  if  any,  models  that  can  even  attempt  to  treat 
~heae  effects  in  a  general  manner  for  snow. 

The  a  parameters  taken  or  derived  from  the  literature  are  shown  in 
figure  1  along  with  the  various  distribution  parameters  used  in  a  corre¬ 
sponding  rain  model.  The  scarcity  of  snow  data  is  evident,  but  the  data 
are  consistent  with  expectation.  Dry  snow,  representing  pure  ice  crys¬ 
tals  and  air,  has  small  MMW  refractive  indices  relative  to  liquid  water, 
and  correspondingly  low  extinction.  As  snow  becomes  wetter,  its  scatter¬ 
ing  cross  section  increases  dramatically— an  effect  seen  as  the  radar 
"bright  band."  At  higher  frequencies,  as  scattering  approaches  the  opti¬ 
cal  limit,  the  cross  section  evidences  a  */R  dependence.  Dry  snow,  with 
its  predominance  of  small  flakes,  has  less  extinction  at  low  frequencies, 
but  as  the  optical  limit  is  approached  the  relationship  of  dry  and  wet 
snow  extinction  will  reverse,  giving  dry  snow  greater  attenuation  than 
wet  snow.  Where  this  crossover  point  occurs  is  not  Indicated  in  this 
data  set.  All  the  curves  pass  through  the  312  GHz  point  for  lack  of  data 
at  these  high  frequencies. 

The  use  of  an  aR*5  relation  for  snow  is  supported  by  this  data,  which 
does  not  display  large  unexpected  variations  from  the  rain  cases.  Note 
that  at  low  frequencies,  below  35  GHz  where  the  Rayleigh  approximation 
applies,  the  curves  for  wet  snow  and  rain  do  not  converge,  due  to  refrac¬ 
tive  index  differences  and  to  the  fact  that  the  extinction  is  expressed 
as  a  function  of  accumulation  rates,  not  mass  density* 

The  b  parameter,  which  determines  the  dependence  of  extinction  on 
accumulation  rate,  is  shown  in  figure  2,  again  with  rain  values  and  the 
inferred  frequency  and  snowflake  behavior. 
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Figure  2*  Snow  parameter  b  with  rain  comparison*  Data  extrapolated 
below  10  CBx 
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_ The  wet  and  nolat  anow  data  are  similar  to  the  Joss  thunderstorm  and 

Marshall-Palmer  behavior;  theae  curves  have  been  roughly  followed  in 
defining  the  wet  and  nolat  anow  b  paraneter  curves .  Dry  snow  behavior  Is 
expected  to  follow  the  Joas  drizzle  distribution  due  to  the  predominance 
of  small  particles  in  both  cases.  However,  two  points  provide  only  the 
crudest  indication  of  the  dry  anow  behavior.  The  large  peak  value  of  -he 
dry  anow  b  paraneter  should  not  cause  alarm,  because  dry  snowrates  are 
generally  quite  low,  there  being  little  or  no  data  for  rates  above  2  ma/h 
in  the  literature. 

Air  temperature  at  the  ground  la  used  to  determine  the  anowtype  for 
an  extinction  calculation.  Dry  anow  la  used  at  or  below  freezing,  moist 
anow  parameters  are  uaed  from  0*C  to  2#C,  and  wet  snow  is  employed  above 
2°C.  These  temperatures  are  rough  estimates,  based  upon  the  assumption 
that  temperature  lanae  rates  are  negative  and  are  subject  to  adjustment 
depending  on the  specific  storm  characteristics. 
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The  frequency  range  covered  by  this  model  is  not  well  defined-  With 
eight  data  points  defining  three  curves  and  with  only  two  measurements 
above  50  GHz,  the  high  frequency  behavior  of  the  model  is  obviously  sub¬ 
ject  to  error.  Por  this  reason,  use  of  the  model  above  100  GHz  should  be 
accompanied  with  caution. 

2.  SNOW  BACKS CATTER 

HHW  anow  backscatter  data  is  severly  limilted  fcy  the  use  of  inferred 
reflectivity  factors  derived  from  measured  snowflake  size  distribu¬ 
tions,  Such  results  are  valid  only  in  the  Rayleigh  regime  below  35 
GHz.  Initially,  a  Rayleigh  model  was  employed  with  restricted  frequency 
coverage.  With  the  publication  of  a  datum  at  95  GHz  resulting  from  the 
SHOW-ONE  measurements,  it  was  decided,  due  to  the  inadequacy  of  the 
existing  model,  to  develop  a  new  model  ualag  ell  the  available  data.  No 
validation  of  this  model  ie  possible  because  all  the  data  were  incorpo¬ 
rated  la  It. 

The  similar  behavior  of  the  rain  and  snow  a  and  b  parameters  led  to 
the  use  of  the  rain  oaekscatter  model  to  provide  the  shape  of  the  snow 
backeciitter  curves.  Thia  procedure  can  be  dangerous,  but  theory  end 
experiment  are  not  yet  able  to  provide  sufficient  knowledge  to  define  an 
independent  model.  At  long  wavelengths,  the  reflectivity  factors  of  Imal 
et  al  (12)  are  uaed  along  with  the  Rayleigh  relation  to  compute  back- 
acatter  cross  section,  n, 
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■nhagfc  k  «  .(a2  -l)Z{m2_  +  2),  m  is  the  complex  index  of  refraction,  X  the 
wavelength,  and  Z  ■  /  N(D)D6  dD  where  N(D)  is  the  snowflake  size  distri¬ 
bution  function.  In  practice,  reflectivity  is  related  to  rain  equivalent 
snowrate  as 

Z  -  aRb  . 


Table  1  lists  the  values  of  the  e  and  b  parameters  from  Iaai  et  al,  where 
their  classification  is  assumed  to  correspond  with  our  dry,  moist,  and 
wet  classes.  - 


TABLE  1.  RAYLEIGH -BACKSCATTER  PARAMETERS 


Snowtype 

a 

b 

1.8 

600 

moist 

1800 

1.8 

wet 

2400 

1.8 

At  shorter  wavelengths,  the  rain  backscatter  curve  was  forced  to  go 
through  the  datum  from  Nem&rich  et  al  (13)  at  95  GHz.  Their  nine  data 
points  were  abstracted  into  a  single  upper  bound  value.  This  was  neces¬ 
sary  because  their  data  did  not  correlate  well  with  the  snow  accumulation 
data,  an  apparent  example  of  snowflake  resonant  scattering  variations. 
The  value  used  was  6.6  x  10"5  m2/m3  for  a  dry  snowrate  of  1  mm/h. 
Examples  of  the  computed  snow  backscatter  cross  sections  for  the  three 
snowtypes  and  a  snowrate  of  5  mm/h  are  shown  in  figure  3.  Comparison 
with  rain  backscatter  indicates  that  for  dry  snow  the  cross  section  is 
much  less  fox  snow;  for  moist  snow,  rain  and  snow  are  comparable;  and  for 
wet  snow,  snow  backscatter  is  much  larger  than  rain. 

In  scaling  the  snowtype  results,  it  was  assumed  that  the  Imai  et  al 
results  would  hold  at  all  frequencies,  thus  giving  moist  snow  three  times 
and  wet  snow  four  times  the  dry  snow  cross  section.  This  nay  not  be 
correct,  but  no  data  are  available  yet  to  improve  the  -viv^ntion. 

3.  MODEL  EVALUATION 

The  SNOW-ONE  test  (3)  of  January  and  February  1981  provided  the 
first  opportunity  to  evaluate  the  snow  extinction  model  using  a  totally 
independent  data  set. 


STAMP  Security  Classification 

here 


STAMP  Security 
Classification 
here 


Insert*  last 


•mme  of  au¬ 
thor!  s)  here 


Start  here 
lor  all 


pages 


after  the 
first 


On  first 
page  type 

title  of _ ^ 

paper  here 

Author  — * 
Affiliation 
City,  State 


first  line 


BROWN 


FREQUENCY 

[«Hz] 


[Figure  3.  Snow  backacatter  cross  section. 

Both  Harry  Diamond  Laboratories  (HDL)  and  the  US  Arny  Ballistics 
Research  Laboratory  (BRL)  carried  out  MKW  propagation  experiments  at 
SNOW-ONE.  The  HDL  Mobile  Measurement  Facility,  operating  at  9$,  140,  and 
225  GHz,  measured  extinction  and  backacatter,  reported  by  Nemarlch!  et  al 
(13). |  The  BRL  experiment  measured  extinction  at  35,  95,  140,  and  217  GHz, 
reported  by  Bauerle  and  Knox  (14).  Due  to  the  small  data  sample  at  35 
GHz,  this  frequency  is  not  included  in  the  evaluation. 

Supporting  meteorological  and  enow  characterization  data  have  been 
provided  by  CREEL  (3)  and  the  US  Army  Atmospheric  Sciences  Laboratory 
(15). 
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The  most  significant  storm  of  the  test,  beginning  at  0400  on 
8  February  and  finishing  at  0500  on  9  February,  mas  selected  for  anal¬ 
ysis.  Snow  accumulation  data  was  available  from  1700  on  tbe  8th  to  0400 
on  the  9th.  At  2200  on  the  8th  windspeede  Increased,  producing  signify 
cant  amounts  of  blowing  snow  and  ending  the  capability  of  the  snow  accu¬ 
mulation  gauges  to  indicate  the  amount  of  airborne  snow*  Therefor#,  only 
the  data  from  1700  to  2200  on  8  February  are  suitable  for  comparison  with 
the  model  predictions.  To  provide  a  sens#  of  the  varletlona  in  extinc¬ 
tion,  the  data  for  the  period  from  1700  on  the  8th  to  0400  on  the  9th 
will  be  presented. 
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’ - The  millimeter  tranamlsalon  data  was  reported  relative  to  the  clear 

path  transmission  determined  either  before  or  after  the  storm.  The  mea¬ 
surements  do  not  correct  for  temperature  and  relative  humidity  variations 
during  the  storm,  but  the  relative  humidity  remained  above  95  percent 
during  this  period,  and  temperature  variations  are  estimated  to  introduce 
at  most  errors  of  ±0.06  dB/km  at  95  GHz,  ±0.15  dB/ka  at  140  GHz,  and  ±0.4 
dB/km  at  220  GHZj — Theae  are  a«all  enough  to  be  neglected. 

The  presence  or  absence  of  fog  during  the  storm  can  Influence  milli¬ 
meter  extinction.  Light  scattering  Instruments  were  used  to  measure 
aerosol  size  distributions,  hut  their  results  are  not  easy  to  inter¬ 
pret.  At  the  succeeding  SNOW-ONB-A  test  held  in  Dec  81  and  Jan  82,  a 
background  aerosol  was  measured  at  all  times,  presumably  due  to  the  wood 
stoves  used  to  heat  many  houses  in  the  region.  At  SNOW-ONE  this  back- 
grouud_aerosQl.  waa.Jiot  anticipated  and  therefore  not  measured.  It  is  not 
-possible- to  separate  unequivocally  the  fog  or  smoke  aerosol  contribution 
to  the  aerosol  distributions  measured  at  SNOW-ONE*  During  the  period 
from  1700  to  2200,  equating  all  the  measured  aerosol  with  fog  provides 
only  a  maximum  of  0.035  gm/m3  of  fog  liquid  water,  and  generally  much 
less.  This  corresponds  to  a  maximum  of  0.17  dB/km  at  95  GHz,  0.24  dB/km 
at  140  GHz,  and  0.29  dB/km  at  220  GHz,  all  less  than  10  percent  of  the 
extinctions  measured  during  that  time.  Apparent  fog  densities  did 
Increase  dramatically  after  2200,  but  this  was  certainly  due  to  the  blow¬ 
ing  snow. 

Bather  thau  compute  the  extinction  from  the  measured  snowrate,  the 
measured  snow  extinction  was  used  In  an  inverse  calculation  to  predict 
the  snow  accumulation  rate.  This  was  done  to  provide  the  moat  concise 
temporal  presentation  of  the  three  frequency  results.  Figures  4  and  5 
\  show  the  ‘‘Inverse'*  model  predictions  compared  to  the  measured  snowrate 
|  for  the  BRL  and  SDL  data,  respectively*  Conversion  from  moist  to  dry 
.  at  1850,  follow iu*  the  temperature  decrease,  has  been  msde.  For  the 
’  u  fro*i  1700  to  2200  the  agreement  is  quite  good  for  an  empirical 
model  with  svrh  a  limited  data  base.  The  peak  at  217  GHz  In  the  BRL  data 
at  21Q0  la  not  repeated  in  the  HDL  data  and  cannot  at  this  time  be 
explained.  Toe  consistency  of  the  model  is  very  good,  snd  It  is  tempting 
to  conclude  that  the  snow  accumulation  measuring  equipment  slightly 
underaaaplea.  Such  a  conclusive  is  unwarranted,  but  does  Indicate  that 
traditional  snowrate  measurements  are  really  not  adequate  for  the 
detailed  characterisation  required  by  this  data. 

The  large  excursions  after  2200  are  due  to  the  blowing  snow,  but  the 
disagreement  between  the  HDL  and  BRL  results  at  225  and  217  GHz  may  be 
due  in  part  to  improper  calibration  of  the  HDL  2 15  GHz  data  becauaa  the 
correct  calibration  was  not  yet  available.  Since  this  anew  model  cannot 
address  the  conditions  after  2200,  discussion  is  academic  at  this  time. 


■try  <jf 

v-r  only  . 


STAMP  Security  Classification 

hero 


4 


tr 


ARMY-  SCIENCE  CONFERENCE  PAPER 


Insert *last 
i.ime  of  au¬ 
thor  (s)  here 

— ► 

Start  here 
for  all 
pages  — ► 

after  the 
tirst 

On  first 
page  type 

title  of _ ^ 

paper  here 

Author  — ► 
Affiliation 
City,  State 


First  line 


STAMP  Security 
Classification 
here 


-BROWN 


.'•TAMP 

.  awngrading/ 
'..'0  lass  i  fi¬ 
at  ion 

.1  format  ion 
i  f  i.rst 
■age  of 
iper  only 


Figure  4.  Snowrate  computed  from  BRL  data  compared  to  measured  snow- 
rate. 
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Figure  S.  Snow  rate  computed  from  HDL  data  compared  to  meaaured  snow¬ 
rate. 
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Theae  results  give  added  confidence  to  the  use  of  the  model  at  fre- 
quencies  up  to  225  GHz,  at  least  for  dry  snow.  Wet  and  moist  snow  accu¬ 
racy  oust  still  be  determined  at  these  high  frequencies.  Given  existing 
field  test  prograas,  it  is  unlikely  that  aeasureaents  above  225  GHz  will 
be  available,  so  that  extension  of  the  model  with  confidence  above  this 
frequency  may  not  be  possible. 


CONCLUSION 


A  new  MMW  snow  aodel  has  been  presented,  based  upon  the  use  of  an 
aRb  relation.  The  data,  used  in  its  develcpaent  were  briefly  reviewed, 
and  the  similarity  of  the  snow  and  ralp  models  was  shown.  The  extinction 
model  was  compared  with  independent  data  from  two  separate  facilities  and 
demonstrated  remarkably  good  agreement.  The  backscatter  aodel  could  not 
be  evaluated  .because  no  Independent  data  exist .  _ _ 
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The  development  of  this  aodel  has  clarified  three  data  deficien¬ 
cies:  lack  of  snow  backscatter  data  from  all  types  of  *nuw,  lack  of 
polarized  propagation  data  at  all  frequencies,  and  lack  of  unpolarized 
data  at  frequencies  above  100  GHz.  Efforts  to  provide  for  the  inadequa¬ 
cies  of  the  data  base  continue,  hut  snow  data  is  difficult  to  obtain,  and 
further  development  of  the  model  will  be  alow. 

A  particular  problem  la  that  a  large  data  base  is  necessary  at  each 
frequency  to  define  the  a  and  b  parameters,  due  to  the  use  of  snowrate  as 
the  fundamental  quantity.  This  also  limits  the  model  to  near  calm  condi¬ 
tions,  since  there  is  no  way  to  estimate  the  effects  of  blowing  snow,  A 
more  sophisticated  model  using  airborne  snow  density  or  the  snowflake 
size  distribution  function  Itself  could  be  developed,  but  the  data  set 
would  be  even  smaller  than  that  employed  here.  Since  these  data  ace  not 
pari  of  the  standard  meteorological  data  base,  it  would  be  moat  difficult 
to  apply  such  models  over  various  geographic  and  diverse  climatological 
areas.  For  these  reasons,  this  empirical  MMW  snow  extinction  aod  back¬ 
scatter  model  is  offered  as  a  solution  to  the  predictive  problem  of  the 
effects  of  snow  on  HMW  propagation. 
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